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ABSTRACT 

We present for the first time metallicity maps generated using data from the Wide Field Spectrograph 
(WiFeS) on the ANU 2.3m of 9 Luminous Infrared Galaxies (LIRGs) and discuss the abundance gradients 
and distribution of metals in these systems. We have carried out optical integral field spectroscopy (IFS) of 
several several LIRGs in various merger phases to investigate the merger process. In a major merger of two 
spiral galaxies with preexisting disk abundance gradients, the changing distribution of metals can be used as 
a tracer of gas flows in the merging system as low metallicity gas is transported from the outskirts of each 
galaxy to their nuclei. We employ this fact to probe merger properties by using the emission lines in our IFS 
data to calculate the gas-phase metallicity in each system. We create abundance maps and subsequently derive 
a metallicity gradient from each map. We compare our measured gradients to merger stage as well as several 
possible tracers of merger progress and observed nuclear abundances. We discuss our work in the context 
of previous abundance gradient observations and compare our results to new galaxy merger models which 
trace metallicity gradient. Our results agree with the observed flattening of metallicity gradients as a merger 
progresses. We compare our results with new theoretical predictions that include chemical enrichment. Our 
data show remarkable agreement with these simulations. 

Subject headings: Galaxies: Abundances, Galaxies: Interactions, Infrared: Galaxies 



1. INTRODUCTION 

Galaxy collisions and mergers represent a key stage in 
the evolution of galaxies in the local universe and beyond. 
Throughout cosmic time gravitational forces have assembled 
ever larger galactic systems from the collisions and mergers 
of smaller fragments. Beginning with the theoretical Toomr^ 
(IT977) sequence of merging galaxies, theoretical modeling of 
massive, merging galaxies now includes detailed physics, al- 
lowing us to make predictions about galaxy evolution over 
cosmic time. The history of chemical enrichment is tied to 
both star formation and the dynamic redistribution of gas 
throughout the lifetime of a gala xy and is dras t ically modi- 
fied by gala xy merger events (iKobavashil l2004t iRupke et alj 
1201 Oat Easy et al. 201 1). 

Theory predicts that major mergers encourage the forma- 
tion of bars in the stellar and gas disks, which induce vig- 
orous gas inflows as the gas looses angular mo mentum to 
the stellar component (Barnes & HemquisBl 19961) . These in- 
flows are thought to be responsible for fueling a massive 
central starburst and feeding AGN and/or quasar activity 
jMihos & Hernquist|[T996l:lBarnes & Hernquisjl 19961) . For a 
spiral galaxy with a preexisting metallicity gradient gas in- 
flow flattens the gradient by diluting the higher abundance gas 
in the central regions with the lower a bundance gas from the 
outer parts of the galaxy (Rupke et al. 2010b a; K ewlev et alj 
l2bl(^) . This flattening is compounded as the spiral arms are 
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Stretched by tidal effects (iTorrev et al.l 1201 lb . Here we in- 
vestigate this process in nearby Luminous Infrared Galaxies 
(LIRGs) 

LIRGs are excellent targets to inform the study of galaxy 
evolution. Although relatively rare, LIRGs compose the 
larger part of the IR luminosity of the universe by z^l 
(ILe Floc'h et al.1 l2005h . As a group LIRGs span merger 
stages from isolated to post-merger and contain massive star- 
bursts, AGN, shocks and a variety of stellar p opulation ages 
(Sanders & Mirabel 1996'; iKewlev et allllOOll: ISanderset al] 
,2003; Yuan e t al. 2010). In most cases the starbursts and 
AGN in LIRGs are driven by ongoing mergers and the pro- 
gression of these mergers is correlated with a rise in the IR 
luminosity, with massive merging LIRGs eventually surpass- 
ing Lju = lO^^i0 to become Ultraluminous Infrared Galax- 
ies (ULIRGs) (e.g. Sanders et al. 1988). Additionally, both 
merger activity and the incidence of U/LIRGs increa se at 
higher redshift (de Ravel et al.ll2009HBundv et al.ll2009l ). It is 
clear that merging, IR-luminous systems played an important 
role during the peak of star formation and in the assembly of 
present-day massive galaxies. This makes LIRGs an ideal tar- 
get to study the mixing and re-distribution of heavy elements 
in mergers to test the pred ictions of merger models. 

IKewlev et alj ( I2006al) found the first observational link 
between merging spiral galaxies and tidally induced gas 
flow in the depressed nuclear metallicities of close merg- 
ing pairs. Other studies of mass-metallicity and luminosity- 
metallicity correlations find that merging syst ems tend to be 
underabundant for their size and brightness ( Lee et al.l 120041: 
IRupke et alj l200l IMson et all 12008) . Theoretical model- 
ing of the metal distribution in merging systems reproduces 
the observed nuclear underabundance and predicts a flatten- 
ing of a preexisting abundance gradient as the merger pro- 
gresses dRupke et ai.li20 1 Oal) . Recent multi-slit spectroscopy 
of HII regions in close-pair spiral galaxies confirmed that gas- 
phase metallicity gradients are indeed flatter in merging sys- 
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terns when compared to a contro l sample of isolated spirals 
(iKewlev et alJl20Tol:lRupke"etan [2010b). 

In this paper we present an integral field spectroscopic (IFS ) 
study of chemical abundances in 9 nearby LIRGs. Our sys- 
tems are a subset of a larger IFS sample of nearby U/LIRGs: 
the Wide Field Spectrograph (WiFeS) Integral Field Unit 
(IFU) Great Observatory All-Sky LIRG Sui-vey (GOALS) 
Sample (WIGS). IFS data provides a wealth of useful infor- 
mation about nearby galaxies, allowing the observer to mea- 
sure several aspects of the gas physics and stars over large 
areas. By using our IFS data to generate maps of emission 
line fluxes, we can similarly create abundance maps and track 
metallicity gradient as a function of merger stage and sample 
later merger stages than previously considered. 

We provide a summary of our sample, the observations and 
the data reduction in Section 2. In Section 3 we discuss the 
analysis of our spectra and our metallicity determinations and 
calibrations. We also present the resulting metallicity maps 
and gradients created from our IFU data in Section 3. We an- 
alyze our observations and compare our measurements with 
previous observations of metallicities in merging systems in 
Section 4. In Section 5 we discuss our observations in the 
context of recent merger models. Section 6 gives our conclu- 
sions. Throughout this paper we adopt the cosmological pa- 
rameters Ho=70.5 km s"^Mpc~\ r^v =0.73, and 17m=0.27 
based on the five-year WMAP resul t s iHinshaw et al.l (l2009h 
and consistent with the lArmus et al.l (120091) summarv of the 
GOALS sample. 

2. SAMPLE, OBSERVATIONS AND DATA REDUCTION 
2.1. Our Sample 

Our targets are drawn from th e Great Observator y All-Sky 
LIRG Survey (GOALS) sample ( lArmus et al.ll2009l) . GOALS 
is a multi-wavelength survey of the brightest 60/xm extra- 
galactic sources in the local universe {log{LjB J Lq) > 11.0) 
with redshifts z < 0.088 and is a comple te subset of the IRAS 
Revised Bright Galaxy Sample (RBGS) dSanders et al.l2003h . 
Objects in GOALS cover the full range of nuclear spectral 
types and interaction stages and are excellent analogs for com- 
parison with high-redshift galaxies. Due to the location of 
WiFeS, our targets comprise a southern sample, with a decli- 
nation limit of about +15°. 

The LIRGs in our sample represent a variety of merger 
stages, nuclear separations, luminosities and physical pro- 
cesses. The systems presented in this paper are those objects 
for which we could measure the extended gas-phase metallic- 
ities with sufficient spatial resolution, as the optical emission 
line gas in these systems is too compact and too extinguished 
to sufficiently measure a metallicity gradient with our ground- 
based data. In addition, we rule out galaxies whose spectra 
are overwhelmingly dominated by non Hll-region emission. 
In our sample th is is primarily widespread radiative shocks 
(iRich et al.ll201 Ih . Systems ruled out consist primarily of the 
post-merger objects in our sample. For our sample, the data is 
insufficient for 10 of the 12 post-merger U/LIRGs in our sam- 
ple: follow-up on larger telescopes with AO and/or Hubble 
Space Telescope would likely allow for a better investigation 
of post-merger targest. 

We analyze a total of 1 1 galaxies in 9 systems, shown in 
Table 1. We classify the m erger stage of eac h system us- 
ing the scheme adopted by lYuan et M (120101) . which is a 
slight modification of the merger stage classification outlined 
by .Veilleux et al.. (.2002.) based on the comparison of obser- 



vations with the simulations of iBames & Hernqui"s3 (119921 
119961) . Isolated systems ('iso') show no signs of interaction 
and have no companion galaxies within 10(3 kpc in projected 
distance. Isolated galaxies may have extremely distant com- 
panions beyond 100 kpc, but they show no morphological in- 
dication of current or past major interaction. Widely separated 
systems ('a') show signs of an ongoing major merger includ- 
ing tidal tails and bridges, but companions are separated by 
a minimum projected distance of 10 kpc. Closely interact- 
ing systems ('b') show more advanced tidal structures and are 
separated by less than 10 kpc, but are not yet coalesced. Fi- 
nally we group all coalesced mergers together ('cde'), includ- 
ing diffuse, compact and post-merger systems. All coalesced 
mergers are characterized by a single nucleus in ground-based 
optical and NIR images, with an increasingly compact core as 
the final stages of the merger take place. Our sample has 2 
isolated systems, 1 in stage 'a', 4 in stage 'b' and 2 in stage 
'cde'. 

2.2. Comparison Samples 

We compare our work to the prev i ous samples 
of IKewlev etalJ ( l20Toh : iRupke et alJ ( l2010bh . iKewley et all 
(l2oTor study a sample of 5 pairs of local luminous spiral 
gala xies from the optically selected galaxy pair sample of 
Barton et al. (2000), each with separations of 15 - 25 kpc. 
Rupke et al.. (.2010b.) expand on this sample of interacting 
systems and also draw from the optically-selected Arp (1966R 
catalog as well the infrared-selected samples of Sa nders et all 
(12003) and Surace et al. (2004) for a total sample of 22 in- 
teracting galaxies in 9 pairs/groups, all classified in stage 'a' 
according to the classification used in this paper ,Rupke et al.l 
( 1201011) also assemble a control sample of local, isolated 
spiral galaxies from the literature with properties similar 
to the spirals in their interacting galaxy sample. Our work 
overlaps with and extends the Rupke et al. sample to later 
merger stages. 

The comparison samples consist primarily of gas- 
rich spirals which are accepted as the progenitors of 
LIRGs (Bai-nes & Hemauist 1996; Mihos & Hernauist 1992 
'Veilleux et al. 2002; lono et al.ii2004tiIshida2004tiNaab et all 
2006; Rupke et al. 2008). The higher average Ljfj of our 
LIRG sample is the result of intense star formation in- 
duced by the merger process-three of the interacting sys- 
tems in the compari son sample are in fact LIRGs themselves 
(IRupke etalJl2010bh . 

Our systems are 1 to 2 magnitudes brighter in M^, which 
generally correlates with stellar mass, as older late-type stars 
trace the majori ty of stellar mass with spectral emission peak- 
ing in the NIR d Aaronsonll 1 9771) . This enhancement is in part 
due to the merger process, as the total Mk of the later-stage, 
combined systems will reflect the total Ma' of the progenitors. 
An additional enhancement in Ma' also reflects a contribution 
from the intense, young starbursts in LIRGs. Both models and 
observations show the possibility of a significant contribution 
to the NIR flux (from 30% to 6 0%) from asymptotic and red 
giant bra nch stars (Mouhcine & Lan9onl2002l;lMaraston et al.l 
120061; ?; iMelboume et al.ii2012i) . Thus we expect the to 
be enhanced of order 1 to 2 magnitudes at approximately sim- 
ilar stellar mass indicating the masses of the systems in our 
sample are roughly consistent with the masses of the com- 
parison sample. This means Ma may not correlate well with 
stellar mass in LIRGs, but the enhancement in Ma should 
still trace the merger process in the same way as Lju. 
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2.3. Observations & Data Reduction 

Our data were taken using WiFeS at the Mount Stromlo and 
Siding Spring Observatory (SSO) 2.3 m telescope. WiFeS 
is a new, dual beam, image-slicin g IFU comrn i ssione d in 
May 2009 an d described in detail bv lDopita etaTI (l2007h and 
iDopita et alj (|2010j). Our data consists of blue and red spectra 
with a resolution of R 3000 & R 7000 and wavelength cov- 
erage of -3500-5800 A & -5500-7000 A respectively. The 
data were taken over 5 separate observing runs in July, August 
and September 2009 and in March and May 2010. 

The d ata were reduced an d flux calibrated using the WiFeS 
pipeline ( iDopita et alj2010h . The pipeline uses IRAF routines 
adapted from the Gemini North's Near-Infrared Integral Field 
Spectrometer (NIFS) data reduction package. A single WiFeS 
observation consists of 25, l"wide, 38"long slit spectra with 
contemporaneous sky spectra if the observation was taken in 
nod and shuffle mode. All of the final combined data cubes in 
our sample consist of at least 2 observations of the same target 
and in some cases a mosaic of two or more pointings. Where 
there were not sufficient observations to remove cosmic rays 
via median-combinati on, cosmic r ay removal was performed 
with the "dcr" routine (lPychll2004h . In this paper we provide 
only a brief summ ary of the d ata reduction, the process is de- 
scribed in detail in iRich et all f2010,2011). 

Individual observations are bias-subtracted using bias 
frames taken as near in time as possible to the observation 
frames to avoid temporal effects. The observations were flat- 
fielded using quarts lamp flats and twilight sky flats. The 
individual spectra are spatially calibrated using a thin wire 
combined with a continuum lamp. CuAr and NeAr arc lamp 
spectra taken throughout each night are used to wavelength 
calibrate each observation. The 25 resulting slitlet-spectra are 
then rectified into blue-arm and red-arm data cubes and sam- 
pled to a common wavelength scale. 

Telluric lines are then removed using observations of B- 
stars or featureless white dwarfs taken at similar airmass the 
same night as the data. Each data cube is then flux calibrated 
using observations of flux standard stars taken on the same 
night. If there were any non-photometric pointings of a target, 
these pointings were scaled and flux calibrated using photo- 
metric data of the same target. We correct for the effects of 
atmospheric dispersion using the WiFeS pipeline. 

Each individual reduced data cube from observations of a 
single galaxy is then combined to form a final data cube sam- 
pled to a common spatial grid for analysis purposes. The 
WiFeS detectors have 0.5"pixels along the slit, so data taken 
in 2010 are binned on-chip by 2 pixels in the spatial direction 
in order to increase signal to noise and produce square spatial 
elements 1" x 1". Observations taken in 2009 are binned post- 
reduction to achieve similar results. Typical seeing achieved 
at SSO during our observations is — 1.5", with some varia- 
tion, on par with the spaxel size for our data cubes. We do 
not perform further binning (e.g. Voronoi tessellation) prior 
to our data analysis. 

Combined data cubes were aligned astrometrically by com- 
paring a pseudo 'r-band' image generated using the red spec- 
trum from each spectral pixel (spaxel) with either DSS R-band 
or HST data where available. Deprojected radii were calcu- 
lated using inclination and position angle data from Hyper- 
Leda (Paturel et al. 2003), with the central spaxel defined us- 
ing the WiFeS 'r-band' images. These data, as well as optical 
diameters are provided in Table 1 . In the case of a few of our 
later-stage mergers we assume an uninchned disk when calcu- 



lating deprojected radii-these have no PA or Inclination given 
in Table 1 . This provides the most conservative estimate of a 
gradient in that system while avoiding ambiguities caused by 
the complex morphologies of late-stage major mergers. 

2.4. Spectral Fitting 

We analyzed every spectrum using an automated fitting rou- 
tine written in IDL, UHSPECFIT, which is based on the code 
created for the work in Zahid et a l. (2011) and is also em- 
ployed b y and described with e xample fits in iRupke et al.l 
(I2010bl) : iRichetalJ JIoToi IMlh . Our routine fits and sub- 
tracts a stellar continuum fr om each spectrum using popu- 
lation synthesis models from ' Gonzalez Delgado et akl (l2005h 
and an IDL routine which fits a linear combination of stel- 
lar templates to a galaxy spectrum using the method of 
Moustakas & Kennicutt ( 2006 a). 

Lines in the resulting emission spectra are fit using a one 
or two-component Gaussian, depending on the goodness of 
fit determined by the routine. All of the emission lines are fit 
simultaneously using the same gaussian component or com- 
ponents. Both continuum and emission lines were fit using 
the MPFIT package, which performs a lea st-squares analysi s 
using the Levenberg-Marquardt algorithm (lMarkwardtl2009l) . 
All of the emission line fluxes used in this paper are subject 
to a minimum signal-to-noise cut of 5. 
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Table 1 


IRAS# 


Other Name 


log(L/ii/L0) 


Ma- 


PAD 


Inclination (°) 


d25 (') 


Merger Stage 


F01053-1746 


IC 1623 A/B 


11.71 


-24.61 


63.4 


34.8 


1.05 


b 


08355-4944 




11.62 


-23.70 






0.30 


cde 


F10038-3338 


ESO 374-IG032 


11.78 








0.76 


cde 


F10257-4339 


NGC 3256 


11.64 


-24.78 


83.2 


48.7 


3.31 


b 


F13373+0105 EAV 


Arp240 


11.62 


-24.96/-25.13 


177.9/85.1 


34.2/62.1 


1.48/1.48 


a 


F17222-5953 


ESO 138-G027 


11.41 


-24.23 


41.7 


51.2 


1.07 


iso 


F18093-5744 N/S 


IC 4687/4689 


11.62 


-24.47/-23.90 


51.2/141.4 


52.7/74.0 


1.15/0.91 


b 


F18341-5732 


IC 4734 


11.35 


-24.48 


102.6 


57.4 


1.32 


iso 


F19115-2124 


ESO 593-IG008 


11.93 


-25.68 






0.55 


b 



Note. — WiFoS GOALS metallicity sample. Names and IR luminosities are taken from the lArmus et al] 1)20091 ) 
summary of the GOALS sample. Mj^ is derived using 2mass data. PA, inclination and d25 for each system are all 
taken from the HyperLeda derived values (Paturcl ct al. 2003) when possibl e, those with suspe ct /unavailable values 
are noted with The merger morphology scheme values (adapted from [Yuan et aO II2010I ) when available), are 

given in the final column. Some of the sources classified as "a" or "b" have multiple pointings of individual galaxies. 



3. ANALYSIS 

3.1. Extinction Correction 

To account for offset in flux between the red and the blue 
data cubes for each system, we compare the balmer decrement 
for Ha, which faUs on the red arm, and H/? on the blue arm 
to the decrement between H/3 and H7 which both fall on the 
blue arm. The average scaling factors were generally within a 
few percent of unity, and were used to scale the red fluxes. 

Due to the low flux and thus limited spatial coverage in H7, 
we did not use H/3 vs. H7 to generate the extinction maps 
for our final analysis. Final extinction maps were created for 
each data cu be using the H a /H/3 decrement and the extinc- 
tion curve of I Cardelli et al.l (IT989i) . These maps were then 
smoothed with a 2-d boxcar size of 3 spaxels to reduce noise 
in the maps. The extinction map for each system is provided 
in the appendix. We then deredden all of the measured emis- 
sion line fluxes for a system using our derived extinction map. 

3.2. Emission-Line Ratios 

We generate emission line ratios from the extinction- 
corrected line fluxes to diagnose the source of excitation for 
each spaxel in the systems we analyze. Emission line ra- 
tio maps are a convenient method of tracing the excitation 
mechanism within galaxies: region of high [N II]/Haor 
[S II]/Ha are generally associated with non-HII region pho- 
toionization. We subsequently classify each spectrum of ev- 
ery individual spaxel as Hll-region like, composite or AGN- 
like (LINER or Seyfert) using the [N II1/Ha . [S II I /Haan d 
[O III]/H/3 ratios and the scheme of 'Kewl ev et alj ( I2006bl) . 
Emission line ratio maps for each system are given in the ap- 
pendix, and are shown without extinction correction applied. 

To avoid contamination from non-HII region photoioniza- 
tion , we reject spaxel s with e mission line ratios that lie above 
the iKauffmann et alj ( 120031) empirical pure star formatiion 
line in the [N II]/Havs. [O III]/H/3 diagram, as well as 
points that lie beyond the theoretical pure star f ormation line 
in the [S II]/Havs. [O III]/H/3 diagram of iKewlev etal] 
(l2006bl) . To account for errors we allow points that fall within 
0.1 dex to the right of th e cutoffs, consistent with the work 
of iRupke et alj fcoiObI) . We do not use the [O I]/Havs. 
[O III]/H/3 diagnostic due to the significantly lower S/N in 
[O I] for most of our systems. 

The spectra that fall into the composite region and AGN 
classes are contaminated primarily by shock excitation in our 
systems. This is apparent in several of the emission line ra- 
tio maps, where off-nuclear regions show enhanced emission 
line ratios. In .Rich et al.. (.201 1.) we showed that these higher 



emission-line ratios are reproduced well by new slow-shock 
models and investigate the star-forming and shocked compo- 
nents separately using kinematic information from our WiFeS 
spectra. In the case of the galaxies IRAS F10257-4339 and 
IRAS F01053-1746, we are able to separate the shocked com- 
ponent from the pure Hll-regio n componen t kinematically us- 
ing this method as outlined in iRich et al.l (1201 lb . For these 
two systems we apply a velocity dispersion cuts of 65 km/s 
and 90 km/s respectively as well as line-ratio classifications 
described above, leaving only pure Hll-region emission line 
ratios for gas-phase metallicity calculation. 

3.3. Abundance Catenations 

We calculate the oxygen abundance in each spaxel using 
several strong-line gas-phase metallicity calibrations. We are 
unable to use the "direct" Tg method due to the low flux of the 
[O III]A4363 auroral line and the higher metallicities found 
in our objects. We instead employ both diagnostics empiri- 
cally calibrated against the Tg method and theoretical meth- 
ods based on the strengths of other measurable emission-line 
fluxes, adopting the procedures ou t hned for each method in 
the summary of iKewlev & EUisonl ( 12008 ^. We measure rela- 
tive abundances within a given diagnostic, avoiding the prob- 
lems associated with absolute ab undance cal culations (e.g. 
iLopez-Sanchez & Esteban 2010; Lopez-Sanch ez et al. 201 11 

We employ the empirical calibration of Pettini & Paggl 
120041 (PP04), which uses the fine ratios [N II]/Haand 
[O III]/H/3 . Using these line ratios has the benefit of avoiding 
the necessity of extinction correction-PP04 developed their 
calibrations for use in higher redshift systems where proper 
measures of the reddening may not always be possible. 

We also calculate metallic ities with the theoretical method 
of IKewlev & Dopital 120021 (KD02), specificaUy we use the 
[N II]/[0 II] line ratio which is insensitive to variations in 
ionization parameter Previous work has shown that the KD02 
diagnostic shows less RMS dispersion than other diagnostics 
for this reason (e.g. Kewlev & Ellison 2008). When making 
comparisons with the iRupke et al. ( 20 lObRand IKewlev et alj 
(I2OIO1) gradient studies, we convert abundances measured 
with the PP04 di agnostics into the KD02 diagnostic using the 
prescriptions of IKewlev & Elhsonl ( 120081) for consistency. 

We also investigated the R23 

([O II]A3727h-[0 III1AA4959, 5007)/ H^ measures of abun- 
dance by employi ng the calibration s of iKobuInickv & KewlevI 
|200l(KK04), and lMcGauglIT99ll (M9I). The KK04 mefliod 
corrects for variations in ionization parameter, by using the 
[O III]/[0 II] line ratio in an iterative fashion to calculate 
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both an ionization parameter and a metallicity. Our results 
using the i?23 method are consistent with PP04 and KD02 
with somewhat increased RMS scatter and are thus excluded 
from further discussion in this paper 

3.4. Metallicity Maps & Gradients 

We use the derived abundances to create maps of the metal- 
licity for each system in our sample. We then measure the 
metallicity gradient as a function of deprojected radius us- 
ing the metallicity values calculated for each spaxel. We cal- 
culate our gradients using an unweighted least-squares linear 
fit, with errors de rived using monte-ca rl o methods, consisten t 
with the work of iKewlev et il] dloToh : iRupke etalJ j2010bl) . 
Only points which are below the Hll-region cut are shown on 
the metallicity maps and gradient plots. This is evident when 
comparing the line-ratio maps to the metallicity maps: regions 
of strong line ratios and/or low surface brightness in the line 
ratio maps are rejected and not shown on the metallicity maps. 
IRAS F18341-5744 and IRAS F13373H-0105 W for instance 
are dominated by strong shock-like emission over large areas 
near the nuclei of those systems. 

We calculate gradients in both dex/kpc and dex/R25, where 
R25 is the effective optical radius, with the latter measure 
meant to account for differences in physical scale from sys- 
tem to system. In general we are not able to e xtend our gra- 
dients t o the same fractional isophotal radius as iKewlej et alj 
(l20Tdh : iRupke et alJ (I2010H ). In some of our systems this is 
due to the spatially restricted, highly extinguished nature of 
the line-emitting gas, coupled in some cases with the limited 
total field of view of our WiFeS mosaics. Our monte-carlo 
error calculation takes this into account, and we also test the 
effect smaller spatial sampling may have on our gradient mea- 
sures by taking systems with spatially well-sampled gradients 
recalculating new gradients at various limited spatial scales. 
The net effect when considering only the inner regions (from 
0.0 to 0.4 R/R25) is the calculation of a slightly steeper gra- 
dient than expected. When other portions (e.g. from 0.2 to 
0.6 R/R25) of the gradient are sampled, the resulting gradi- 
ent measured is within a few percent of the expected gradient. 
These results indicate that our conclusions regarding gradient 
flattening remain relatively unaffected by spatial sampling. 

3.5. Calibration Discrepancies 

We calculate metallicities for each spaxel using different 
methods in order to investigate any discrepancies in metallic- 
ity measurement introduced by errors in our extinction cor- 
rection, flux measures and from ionization parameter varia- 
tion. Due to the wide separation in wavelength of the diag- 
nostic emission lines the KD02 calibration is most sensitive 
to errors in extinction correction while the PP04 diagnostic 
is essentially immune to reddening effects due to the close 
proximity of the line ratios employed. On the other hand, the 
PP04 calibration should be the most sensitive to variations in 
ionization parameter, while the KD02 diagnostic is essentially 
insensitive to ionization parameter 

In the case of very dusty LIRGs, the PP04 diagnostic proves 
the most useful as it avoids using the very blue and thus very 
extinguished [O II] line. Indeed the PP04 diagnostic provides 
the least RMS scatter from a straight line fit in our gradient 
measures, owing most likely to errors in our extinction cor- 
rection which increase the scatter in O/H values calculated 
using KD02. The inability of the PP04 diagnostic to account 
for variations in ionization parameter, however, is a danger if 



the parameter varies in any systematic way, especially with 
radius. We find a correlation between ionization parameter, 
as measured using the KK04 method, and the difference be- 
tween metallicities measured using the PP04 diagnostic and 
the KD02 or KK04 diagnostic. Fortunately, the effect of vary- 
ing ionization parameter is minor and does not affect the con- 
clusions in this paper Overall, the discrepancies between the 
metallicity calibrations in our analysis are consist ent with the 
comprehensive study of iKewlev & EllisonI (120081) . 

It is also worthwhile to note that our gradients are de- 
rived using integrated spectra containing all of the light 
within each spaxel, including any diffuse emission not as- 
sociated with HII regions. This is in contra st to the tar- 
geted H ll-region abundance gradients used in iKewlev et alj 
( 120 1 01) : IRupkeetalJ (l20I0bh . Although we carefully re- 
move non-photoionized spectra, we might still expect a dif- 
fuse contribution to line-ratios used in gradient calculation 
(A lonso-H errero et al. 2009, 2010). Studies have shown, 
however, that abundances inferred from an integrated spec- 
trum are consistent with Hll-region abundances, regardless 
of the abundance calibration used dMo ustakas & Kennicuji 
'2006bt iRosales-Ortega et al.l 1201 II) . pRosales-Ortega et all 
(201 1) in particular present highly resolved IFU spectroscopy 
of nearby NGC 628 and show that the abundance gradient de- 
rived from their total integrated spectra is very consistent with 
the gradient derived with spectra isolated from Hll-regions. 

4. ABUNDANCE GRADIENTS 

IKewlev et al] |2010') and Rupkeetin ( 12010 b') found a flat- 
tening in the metallicity gradients of widely separated pairs of 
galaxies. Their sample represents systems caught in the earli- 
est stages of merging between first and second pericenter As 
the galaxies continue to merge, gas should continue to flow 
towards the individual galaxy nuclei, fueling the ongoing nu- 
clear starburst and quenching star formation in the tidal rem- 
nants. The nuclear metallicity would first see a depression due 
to the infalling gas, followed by a subsequent enrichment and 
further depression as the merger progresses to coalescence. 
The metallicity gradient is continuously flattened by the com- 
bined effects of infalling low-metal gas and tidal stretching of 
the spiral arms and outer portions of each system. In this sec- 
tion we show that our results indicate the trend of flattening 
metallicity gradient may extend to later merger stages. 

4.1. Merger Stage and Gradient 

In the context of the lYuan et aP (|2010') merger scheme, the 
control and i nteracti ng samples of 'Kewley et al. (2010) and 
iRupke et alJ (120 1 Obi) are respectively 'isolated' and widely- 
separated 'a' stage systems, while the majority of the galaxies 
in our sample are close-pairs in the 'b' stage, further along 
in the merger process. Our sample also has 3 late-stage coa- 
lesced systems ('cde'). It is not clear that all of the stage 'b' 
LIRGs in our sample will becom e true ULIRGs given their 
morphology and Lju. I Yuan et aP (12010) suggest that systems 
with mass ratios closer to unity are more likely to become true 
ULIRGs during the later stages of the merger process. This 
could account for systems such as IRAS F2 1330-3846, which 
has a comparatively low total Ljji despite the fact that it is 
a more advanced merger The gas fraction and dust mass of 
each system prior to the merger also strongly affect the total 
IR luminosity. 

We plot metallicity gradient as a function of merger stage 
in Figures 1 and 2. Although we do not have a large control 
sample of isolated LIRGs, the metallicity gradients of our 2 
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Figure 1. Metallicity gradient as a function of merger stage Oiisolated, liwidely separated pairs ('a'), 2:close pairs ('b'), 3 : coalesced systems (' cde'). The 
isolated control sample of Rupke et al. (2()10b) is represented by the blue circles, the interacting sample of IKewlev et aljf2Q10D : IRupke et an)201Qbl) by yellow 
and our data is plotted as red circles. Simple linear regressions are overplotted for each pair of quantities. 
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Figure 2. Same as figure 1, but plotting dex/kpc instead of dex/optical radius (R25). 
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Table 2 



Variables Correlation Coefficient Probability (%) 



Mk,PP04, R25 -0.49 0.2 

Ma',KD02, R25 -0.48 0.3 

Ljjj, PP04, R25 0.56 0.02 

Ljfl, KD02, R25 0.61 0.005 

Ma , PP04,Rfcpc -0.60 0.008 

Ma, KD02,Rfepc -0.59 0.01 

LiR, PP04, Rkpc 0.65 0.001 

LiH, KD02, Rkpc 0.66 0.0006 



Note. — Spearman rank correlation coefficient. Two- 
sided probability of finding correlation coefficient by 
chance. 

isolated LIRG systems fall within the range of gradi ents of 
the con trol sample of local, isolated spirals from Rup ke et alJ 
(12010b) . Figure 1 shows that metallicity gradient flattens 
through each merger stage, with the three coalesced systems 
showing the flattest gradients. The average gradient in our 
progressed mergers is similar to the average gradient mea- 
sured i n the widely separate d pairs sample of iKewley et al.l 
(l2010l) : lRupke et alJ (l2010bl) . Figure 2 shows that our results 
are unchanged when when metallicity gradient is plotted as a 
function of optical radius and when it is plotted as a function 
of absolute deprojected radius in kpc. 

Interestingly, the galaxy IRAS F13373-I-0015 may be show- 
ing gradient-flattening effects due to both gas infall as well 
as the stretching of its spiral arms due to tidal effects from 
the merger as evidenced by an apparent kink in the metal- 
licity gradient. Unfortunately, the nuclear region is strongly 
contaminated by shock-like emission, so we excluded gas- 
phase metallicity measurements from our plot for these radii. 
If we ignore this contamination and calculate an abundance 
for the center-most region, the measured nuclear [O/H] does 
appear depressed with respect to the observed and calculated 
gradient, which would create a gradient that shows a slightly 
steeper portion between two flatter regions. 

4.2. Luminosity and Gradient 

We also plot both Lm and Mks, quantities which should 
be correlated with merger stage, as a function of metallicity 
gradient in Figures 1 and 2. L//j traces star formation rate 
and is correlated with merger stage, while Mks is unaffected 
by the high extinctions encountered in LIRGs and is a func- 
tion of stellar mass (though see s ec. 2.2). Values for L//? 
are taken from Armus et al. (2009), which adjusts the original 
RBGS CSanders et al. 2003) values to the 5-yeai- WMAP cos- 
mological va lues. Mks is derived using the distances from 
lArmus et al.l (l2009i) and the 2mass 20 mag/arcsec'^ isop hotal 
k- band magnitudes (c onsistent with lRupke et alj (1201 Obi) ). 

iRupke et al] ( 12010 b) concluded there was no strong evi- 
dence of a correlation betwee n gradient and either hjn or 
Mks, though they note that IZaritsky etaH (119941) pointed 
to possible correlations when considering dex/ kpc. When 
our observations are combined with the work of iRupke et al.l 
(1201 Ob) , the case for a correlation between gradient and both 
Lju and Mks correlations is strengthened. The Spearman 
rank correlation coefficients indicate a moderately strong cor- 
relation in all cases, though Lj n is more strongly correlated 
than y^Ks- 

Gradients calculated in dex/kpc are more tightly correlated 
than those measured in dex/R25 (Table 2), consistent with 
IZaritsky et al.l (119941) . Increased scatter when considering op- 
tical radius may be due in part to the increasingly disturbed 




Figure 3. This figure shows the scaled Luminosity & Nuclear abundance 
(L — Z) relation as provided in Fig. 11 of Rupke et al. (2010b). Our work 
extends the trend towards lowered nuclear metallicity for higher galactic lu- 
minosity. This is consistent with the depressed nuclear metalHcity which is 
expected to be seen in merging systems. The data of Rupke et al. 12 0103) is 
plotted in blue (isolated) and yellow (wide pairs) while our LIRGs are plot- 
ted in red. Nuclear metallicity is measured by extrapolating th e gradient to 
Q.IR2 5 and scaling to match the approximate L-Z relation of ISalzer et al J 
(2003) per the method of Rupke et al. (2010b). 

morphologies of merging systems. As a merger progresses 
R25 represents a less consistently comparable value as tidal 
features and the coalescence of both galaxies alter this mea- 
sure such that it is no longer consistent with R25 measured in 
pre-merger spiral galaxies. 

4.3. Nuclear Abundance and Luminosity 

To further compare our sample to previous work, we con- 
sider the nuclear abunda nces from our s ample. In order to 
directly compare to Ru pke et al.l (1201 Obi) we also consider a 
Luminosity-Metallicity {L — Z) relation, rather than a Mass- 
Metallicity (M — Z) relation, as plo tted in fig. 3 . We con- 
sider the the nea r-IR L-Z rela t ion of ISalzer et al.l (120051) for 
consistency with Ruoke et al. (2010b). We adopt the same 
method as Rupke et al. (2010b) for nuclear abundance de- 
termination: we assign an abundance by extrapolating the 
calculated gradients to a fiducial radius of 0.1 R25. We 
also adopt the same offset in absolute metallicity with re- 
spect to the ISalzer et al.l ( 12005b relation. Our nuclear metal - 
licities are indeed lower than the expected L-Z relation, 
again extending previous work on merging systems to higher 
NIR luminosities. This result is consistent with previous 
observations parrs feewlev et al. 2006a; Rupke et al. 200^ 
lEUison et al."2008'; Michel-Da nsac et aL,2008.: ,Peeples et at) 
2009; Sol Alonso et al. 2010) as well as the merging galaxy 
models we consider in our discussion of gradient evolution 
with merger stage and properties. 

5. DISCUSSION 

Our observations are a direct measure of the redistribu- 
tion of gas within strongly interacting systems. The observa- 
tional reality of metallicity gradients in non-int e racting galax- 
ies is well established (Zaritskv et al.l [T994t Ivan Zee et al.l 
119981) . Observations of a change in an established abun- 
dance gradient and systematically lower nuclear metalllic- 
ities are clear indicators of gas flows and have been seen 
not only in strong interactions but also in bar structures 
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Figure 4. Merger stage vs metallicity gradient as measured from our theo- 
retical merger models. The gradients measured from the merger models in 
dex/kpc show values and flattening consistent with observations. The aver- 
age value for the models generated is given by the large black point, the error 
bars represent the spread in model values. Our gradient data from figure 2 
is oveiplotted for reference, data from Rupke et al. (2010b) is overplotted in 
blue points (isolated) and yellow points (wide pairs) and our LIRG data is 
overplotted with red points. 



(e.g. IViia-Costas & Edmuiidil [1991 IDutil & Rovl [19991: 
iMartin & R oy 1995). While iFriedli et al.1 (ll99 4Vdemon- 
strated the redistribution of gas traced by metallicity in barred 
galaxies, only recently have numerical simulations been used 
to address the detailed physical motions of gas and metals 
wit hin interacting galaxies. 

'RuDke et al.' feolOa), 'MontuorietaD (1201 Oh and 
U)rrey et al. (2011) presented the first attempts at using hy- 
drodynamical simulations to address the effects of major 
galaxy mergers on nuclear abundances and metallicity gra- 
dients. The three separate investigations all involve mergers 
of massive spiral disk galaxies with 1:1 mass ratios. As noted 
by l^pke et al. (2010a), these conditions are typical for the 
ULIRG formation scenario and are thus a useful direct com- 
parison. 

The models of iRupke et alj (l2010al) represent the simplest 
set of simulations analyzed for metallicity gradient evolu- 
tion both in merger parameter space and in model sophistica- 
tion. Their models use the methodology described in BarnesI 
([2004), but without any ongoing star formation and thus with 
no chemical enrichment. This approach clearly captures the 
effect that low meta llicity gas inflows wi ll have on an evolving 
nuclear metallicity. iRupke et~an ( l2010ab conclude that the nu- 
clear metallicity can be substantially depressed by gas inflows 
and that metallicity gradient can be dramatically flattened fol- 
lowing tidal tail formation. They find, not surprisingly, that 
the depression in the nuclear metallicity is correlated with the 
gas mass that has rn i grated to the nuclear region. 

iMontuori et al.l (120101) empl oy simulations from th e 
GalMer database, described in (IChilingarian et al.l 120101) . 
Their simulations include chemical enrichment due to ongo- 
ing star formatio n. A s a resu lt of including chemical enrich- 
ment, Montuori et al.l (1201 Oh find that the nuclear metallic- 
ity changes non-monotonically as a function of the merger 
as chemical enrichment partially offsets the effect of low- 
metallicity gas inflows. 

Torrey et al.. (.201 1.) have studied merger induced nuclear 



metallicity evolution while considering star formation and 
chemical enrichment over a relative ly wide rang e of me rger 
parameter space. The analysis in iTorrev et al.l (1201 1 ) ex- 
phcitly differentiated the effects of initial metallicity gradi- 
ent re-distribution from ongoing metal enrichment on the nu- 
clear metallicity evolution - allowing for a basic template for 
merger nuc lear metallici t y evol ution to be presented. The sim- 
ulations of Torrey et al.l (1201 ll) included systems with varied 
initial mass and several merger orientations - allowing for a 
more complete comparison to observations. 

5.1. New Models 

The interpretation of the data presented in this paper can 
benefit by considering a realistic theoretical model for the 
evolution of interacting galaxies. In principle, one could use 
numerical models to directly simulate the merging systems 
that have been presented in this paper However, this would 
require accurately setting a large number of parameters for 
each merging system (e.g., the initial stellar mass of each 
galaxy, the initial gas fraction, the initial morphology, the 
merging orbital orientation, etc.). Fortunately there are certain 
aspects of major mergers, (e.g. the formation of tidal features) 
which are generically associated with strong tidal encounters 
and do not sensitively depend on the detailed characteristics 
(e.g. the gas fraction) of the merging galaxies. 

In the context of metallicity gradient evolution, we can 
learn a lot about the expected evolution of a population of 
merging galaxies by exploring a limited number of carefully 
selected merger simulations. Our goal is to present a realis- 
tic theoretical model to augment the interpretation of our ob- 
servational data. To achieve this, we build a limited merger 
simulation suite that we can explore and understand in detail. 
Although this approach will leave certain areas of merger pa- 
rameters space unexplored, we still expect the characteristics 
and tends present in the simulations to match that of the obser- 
vations which allows us to physically probe the driving forces 
behind the observed evolving metallicity gradients presented 
in this paper 

We compare our observations to a suite of galaxy merger 
simulations (Torrey et al. 2011) carried out using GADGET- 
2 (Springel 2005). Our simulated merger suite consists of 32 
simulations which are a chieved by u s ing 16 merger orienta- 
tions (orientations a-p in lTorrev et al.l (1201 lb ) with 2 progeni- 
tor galaxies. The two progenitor galaxies for our simulations 
both have equal total system mass, but with slig htly differ- 
ent initial disk profiles. Both systems include a 'Hemquisa 
^990) dark matter halo (Mh 5.1 x 1G"A/0, a = 22 kpc) 
with rotationally supported exponential gas and stellar disks 
(^rf,i = 3.2 kpc, Rd,2 = 4.6 kpc). The initial gas fraction of 
the system is 25%, which will decrease with time due to star 
formation. 

A substantial discussion of the kind of models used for 
comparison in this paper as well as further analysis of the- 
ory de rived from those models can be found in T orrey et al.l 
(1201 ll) . Here, we briefly describe our methods for tracking 
the nuclear metallicity and metallicity gradient, which serves 
as the key component of our simulation analysis. Our simu- 
lations track a single metallicity scalar value that is allowed 
to increase as a result of chemical enrichment from star for- 
mation. To compare our simulations to the observations, we 
focus only on the Oxygen abundances, which we assume to 
be 30% of the total metal mass. Because most Oxygen is 
produced in type 11 supemovae which return their mass and 
metals to the ISM over relatively short (i.e. ~ 10^ year) 
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Figure 5. Metallicity gradient co mpared to Lj j ;, nucle ar separation and nuclear metallicity. The color contour plots are generated using measurements from 
the merger models while data from'Rupke et al.' CT 10b( ) is overplotted in blue points (isolated) and yellow points (wide pairs) and our LIRG data is overplotted 
with red points. These plots show good agreement between models and observations while providing infoiTnation about how these quantities roughly track with 
merger progress. L/j{ derived from the simulations is calculated based on the instantaneous star formation rate using the formulae of .Kennicutt et al., t2Q09.) . 

metallicity for all gas particles within a central Ikpc sphere of 
the galaxies center (defined by the potential minimum). Or, 
to determine the metallicity gradient, we calculate the aver- 
age metallicity within a series of concentric spherical shells 
centered on the galaxies potential minimum, and find the best 
linear fit of the resulting radius-metallicity data points. 

For an instructive comparison to the observed systems pre- 
sented in this paper, we track the nuclear separation, nuclear 
metallicity, metallicity gradient, and IR Luminosity during the 
merger Using the data from all simulations, we can then form 
an expected evolution and distribution of galaxy properties 
during the merger sequence. The results are plotted in figures 
4, 5 & 6. 
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Figure 6. L/^ vs. projected separation measured from the merger models 
as plotted in fig. 5 with our data overplotted as red points and the widely 
separated systems of RuDke et al. (201()b) oveiplotted in yellow. Despite the 
agreement between observations and models this plot shows projected sepa- 
ration as a less helpful measure of merger progress. 

timescales, we assume instantaneous metal enrichment which 
we practically incorporate into our simulation by setting the 
rate of metal production proportional to the local star forma- 
tion rate. 

We impose metallicity properties on our galaxies at the be- 
ginning of the simulation such that they conform with obser- 
vations of galaxies in the local universe. Namely, we initialize 
the nuclear metallicity of our galaxies onto the sloan digital 
sky s urvey derived mass-metallicity relation (Tremonti et al] 
120041) with the radial metallicity gradient following a decay- 
ing exponential profile with a characteristic leng th propor- 
tional to the disk scale length (iZaritskv et al.lll994h . 

Using this model, the metallicity in our galaxy is deter- 
mined by the re-distribution of the initial metallicity gradient 
due to bulk movements of gas induced by the strong gravi- 
tational interaction and chemical enrichment associated with 
starburst driven star formation. At any time in the simula- 
tion, the metallicity can be calculated for a region of space 
by finding the mass weighted average metallicity of all gas 
particles in the specified volume. For example, we calculate 
the nuclear metallicity by finding the mass weighted average 



5.2. Merger Progress Tracers 

In Figure 4 the merger model metallicity gradient is plotted 
as a function of merger stage where the simulated mergers are 
classified using the same formalism as our observational anal- 
ysis. The solid black circles show the mean simulated gradi- 
ent value of all merging systems at a given merger stage with 
the error bars representing the l-cr standard deviations derived 
from the spread in gradient values for the simulated galax- 
ies in that bin. By inspection, we find that both the models 
and observations show a clear trend toward shallow metallic- 
ity gradients with increasing merger stage. In the models we 
find a strong change in the gradient that occurs immediately 
after first pericentric passage, when the formation of tidal fea- 
tures stretches the initial metallicity gradient and gas inflows 
begin to flood the nuclear region with low metallicity gas. The 
gradient flattening continues as the merger progresses and by 
merger stage cde (i.e. coalesced systems), both the models 
and observed systems have metallicity gradients consistent 
with being flat. 

Figure 5 shows the evolution of the metallicity gradient as a 
function of IR Luminosity, projected separation, and nuclear 
metallicity. L/jj derived from the simulations is calculated 
based on the instantaneous star formation rate using the for- 
mulae of Kennicutt et al. (2009). Each panel shows the dis- 
tribution of simulated galaxies — as the background colored 
contours — along with the observational data. The colored 
contours represent the local density of simulated data points 
(where each snapshot from each merger counts equally as one 
data point) with the solid black lines denoting regions that 
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contain 50, 70, 80, and 90% of the simulated data points. 
As discussed above, we find that the gradient becomes in- 
creasingly shallow as the merger progresses. We also find 
the IR Luminosity increases and the nuclear metallicity be- 
comes slightly depressed. This is not unexpected, as the same 
fundamental process is driving the evolution of all of these 
quantities. Strong tidal encounters that are associated with 
the merger sequence will form tidal features that will stretch 
or flatten the metallicity gradient. The same tidal encounters 
are responsible for driving gas into the central region, which 
lowers the central metallicity and enhances starburst activity. 
Therefore, we expect that the association of flat metallicity 
gradients with enhanced IR luminosity and depressed nuclear 
metallicity is a robust byproduct of the merger process. 

Figures 5 and 6 show the metallicity gradient and nu- 
clear metallicity against observationally projected separation. 
Again we find that the observations and merger models re- 
main in general agreement, but the projected separation is not 
as useful for tracing merger progress. We find the intermedi- 
ate and late stage mergers are closely clustered in projected 
separation. This is a result of the non-monotonic evolution 
of the nuclear separation with merger stage. Statistically an- 
alyzing large surveys should indeed show enhancements in 
the star formation rate, depressions in the nuclear metallicity, 
and a general flattening of the metallicity gradient that evolves 
with projected separation. However, separating intermediate 
and late stage mergers cannot effectively be carried out using 
only the projected separation information. 

5.3. Further Modeling Considerations 

These conclusions will remain fundamentally unchanged if 
we applied a more extensive merger simulation suite based on 
the physical justification that the formation of tidal features 
and driving of gas to the galactic central region - the physi- 
cal engines behind the observed trends presented in this paper 
- are events that are generically associated with the merger 
process. The formation of tidal t ails in particular is a wel l 
known result of a tidal encounter (iToomre & Toomrelll972b . 
which can depend on the merger orbital parameters (e.g. 
iD'Onghia et al. 2010). If the tidal encounter is sufficiently 
weak because the perturbing galaxy is not massive enough or 
too distant, then we will not expect strong tidal features to 
form and we would not expect to see a dramatic flattening 
of the metallicity gradient. However, in the regime of strong 
tidal encounters, we expect tidal tails to form and the metal- 
licity gradient to flatten without any sensitive dependencies 
on, e.g., varied gas fractions or the presence of a bulge. 

The same arguments apply when considering the driving of 
gas into a galaxy's central region. The tidal encounter sets 
up mis-aligned bars in the stellar and gaseous components al- 
lowing the gas to lose angular momentum and fall into the 
nuclear region (Barnes & Hernquist 1996). Additional con- 
siderations do exist, such as the stabilizing effect of stel- 
lar bulges against axis-symmetric perturbatio ns during distant 
tidal encounters (iMihos & HernquisB 119961) or the reduced 
torquing efficiency of very gas rich systems (iHopkins et al.l 
[2009). However, neither of these effects will prevent a 
strong gravitational interaction from driving gas inflows lead- 
ing into fin al coalescence (see, e.g., Mihos & Hernquist 199S 
IHopkins et al. 2009, for detailed justifications). 

These comparisons provide a useful measure of the accu- 
racy of merger models and strengthen the use of such models 
to make predictions about various aspects of galaxy mergers 
and their effect on galaxy evolution. A more thorough discus- 



sion of the kind of models used for comparison in this paper as 
well as further analysis of theory d erived from those models 
can be found in lTorrev et al.l (1201 ll) . 

When our new observations are combined with previous 
work, there is also evidence of a correlation between metal- 
licity gradient and both Ljj^ and Mk , quantities that loosely 
trace merger progress. Assigning a precise merger stage is 
non-trivial. Determining a rough estimate of the time elapsed 
since first pericentric passage for instance is not possible 
when considering only the morphology and projected sepa- 
ration of a merging pair. One solution could be t he use of 
software like Identikit (iBames & Hib bard 2009; B amesl2011l) 
to constrain merger parameters using observed velocity infor- 
mation. A more reliable merger stage could be assigned once 
an accurate model is obtained, allowing for a more detailed 
study of the change in various physical quantities as a merger 
progresses. 

5.4. Presence of Shocks 

There are other factors which hamper gradient measure- 
ments as mergers progress: Yuanetal. (2010) note a 
marked increase in AGN/LINER influence on optical spec- 
tra as a f unction of L j ^ and merger stage, a n effe c t also 
seen by lArmusetall (119901) : IVeilleux et ai] (119951) : iGotol 
(120051) : lEllison et al.l (120111) . This effect, coupled with the 
increasingly compact region of star formation as the remain- 
ing gas in a merger is concentrated entirely within smaller 
nuclear regions makes measurement of spatially resolved 
gas-phase metallicities difficult and prone to contamination 
from an AGN. In addition, increased contamination from 
radiative shocks would be expected in the most IR lumi- 
nous systems, consistent with increasing gas inflows and the 
obser ved increase in sup ergalactic winds as a function of 
LiR (iRupkeet all 120051) . In our overall IFU sample of 27 
U/LlRGs, ~2/3 of the systems show a contribution from ex- 
tended shock excitation, fractionally increasing as the merger 
progresses. IFU observations of sufficiently high spatial and 
spectral resolution should be able to overcome these issues as 
long as there is of measurable star formation. 

6. SUMMARY AND CONCLUSIONS 

We have measured the gas-phase oxygen abundances in 
9 nearby LIRG systems using integral field spectroscopy. 
We measured emission line fluxes for all of our IFS data 
cubes and used this information to create maps of excitation 
and metallicity measured with various strong-line metallic- 
ity diagnostics. We calculated radial abundance gradients us- 
ing our metallicity maps coupled with radii and deprojected 
radii calculated for each system. Our results agree with re- 
cent observations o f gradient flattening in interacting pairs 
dKewlev et alj|2010l: iRupke et alj|2010bf) and indicate some 
evidence of flattening at later stages of the merger process, 
though further observations are needed to confirm this result. 

Our results are consistent with recent nume rical models 
of metalhcity and gas flows i n major mergers (iRupke et al.l 
1201 Oat lMonfiioriet aLll2010l: iTorrev et al.l 1201 lb . We com- 
pare our observations with a new set of merger models used to 
track various quantities including metallicity gradient, nuclear 
metallicity, Ljji and separation. These new models track gra- 
dients and enrichment carefully throughout the merger pro- 
cess, including the effects of ongoing chemical enrichement 
for the first time. Our observations agree remarkably well 
with the models, adding observational support for theoretical 
predictions about gas flows in numerical merger models. Our 
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comparison between model and observation also indicate the 
efficacy of metallicity gradient and Lju in tracking merger 
progress, while showing nuclear separation as less useful. 

The latest stages of a major merger event prove the most 
difficult period to analyze and interpret using our observa- 
tional methods and data. Investigating the gradient evolution 
using numerical simulations of post-merger systems requires 
fully cosmological simulations where accretion of prestine 
gas at late times can help reestablish an abundance gradient 
(iTorrey et al.li201 1). Investigation of this phenomenon also 
requires a finer spatial scale and deeper observations than we 
have achieved with our current data set. An increasing con- 
tribution from radiative shocks and AGN to the emission line 
spectra of late stage mergers may inhibit these measurements, 
necessitating non emission-line abundance diagnostics. 
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A(dex/R25) PP04 
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08355- 
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0.015 ± 0.018 
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F10038- 
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0.007 ± 0.004 


0.002 ± 0.005 
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9.015 


F17222- 


-5953 


-0.701 ±0.017 


-0.622 ±0.027 


-0.047 ± 0.002 


-0.041 ± 0.002 


9.133 
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F 18093 - 


-5744 N 


-0.259 ±0.023 
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-0.019 ±0.002 


-0.019 ±0.005 
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-5744 S 
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9.045 
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8.876 


Note. 


— Derived metallicity gradient slopes and intercepts in dex/kpc 


and dex/R,25 for 


our sample. 





APPENDIX 
NOTES ON INDIVIDUAL OBJECTS 
Isolated Systems 

There are 2 isolated systems in our sample. These systems have properties most consistent with the lRupke et al.l (1201 Obi) control 
sample, with well established metallicity gradients not yet disturbed by interactions. Both systems are lower luminosity LIRGs, 
withL/H - 1O"-4L0. 

IRAS F17222-5953 (ESO 138-G027). — This system is more akin to a typical, non-interacting, strongly starbursting spiral 
galaxy in our sample. It is in the vicinity of a few other galaxies, including the similarly bright ESO 138-G026, but even the 
nearest galaxy is at a projected distance of over 100 kpc and IRAS F17222-5953 is not yet interacting with any of these systems. 
Our IFU data applied to the [N II1/ Hav [O III1/H/3BPT diagram show a clean curve following the shape of the SDSS sequence 
of local star forming galaxies (e.g. 'Kewley et al.' ('2006b')), but with a slight apparent shift in total [N II]/Ha . We interpret this 
shift as an overall Nitrogen enhancement (Perez-Montero & Contini 2009). 

IRAS F18341-5732 (IC 4734) Like IRAS F17222-5953, our second isolated galaxy is in the vicinity of a few other luminous 
galaxies, but is not yet undergoing any interactions. Our nuclear spectra are dominated by a LINER combined with an aging 
stellar population. The strongest sights of star formation are where the bar in this galaxy meet the spira l arms eviden ced by 
the two strong clumps of Hll-region like points seen in the line-ratio and metallicity. Ha imaging by Dopita et alj (l2002h shows 
further knots of star formation along the spiral arms and our nuclear spectra also show signs of an aging stellar population in the 
nucleus ofIC 4734. 

Widely Separated Systems 

The interacting sample of i Rupke et~an (HoTO^) consists of widely separated systems. Our sample has only a single pair of 
galaxies at this stage of interaction, the interacting system Arp 240. 

IRAS F13373+0105 W (NGC 5257) Although NGC 5257 and its equal mass partner NGC 5258 are still widely separated 
and retain much of their structure, they exhibit interaction features including tidal tails and a bridge between the two galaxies. 
The spiral arms show regions of strong star formation with some signatures of post-starburst populations away from the nucleus. 
The nuclear regions of NGC 5257 are dominated by an older stellar population as the intense nuclear starburst associated with 
the later stages of major mergers has not yet begun. The extinction is higher in the nuclear regions and the measurable line- 
ratios place the nucleus in the composite region of the standard diagnostic diagrams, indicating possible LINER activity-previous 
nuclear obs ervations an d integrated spectrophotometry are consiste nt wit h the higher nuclear extinction and overall line ratios we 
observe (IVeilleux et al]il995l:TKewley et alMOOlUMoustakas & Kennicutl2006al) . 

IRAS Fl 3373+0105 E (NGC 5258) NGC 5258 is a near-twin to NGC 5257 in mass and luminosity. Our spectra indicate higher 
extinction and a flatter gradient, though the extent of measurable Hll-region metallicities is smaller than in NGC 5257. Again 
similar to NGC 5257 the nuclear region of NGC 5258 is dominated by older stars and very little line-emission and there is some 
evidence of younger post-starburst populations away from the nucleus. The strongest line emission is associated with the knot of 
star formation to the southwest of the nuclear region. 

Closely Interacting Systems 

Our sample has several closely interacting systems, we calculate gradients for 5 of them in this paper These systems span a 
range of L/i? from typical LIRG luminosities to nearly ULIRG. IRAS F01053-1746, IRAS F10257-4339 and IRAS F191 15-2124 
are the most indicative of this class, while IRAS Fl 8093-5744 N/S is a triple system with two widely separated spirals and a 
third, less massive closely interacting component. Although IRAS F01053-1746 and IRAS F10257-4339 are both composed of 
two individual galaxies, they have progressed along the merger sequence to the point where we are unable to effectively assign 
gas metallicities to either of the galaxies. We thus treat both IRAS F01053-1746 and IRAS F10257-4339 as individual systems 
with the brightest optical nucleus defined as the zero point for a metallicity gradient. 
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IRAS F01053-1746 (IC 1623) This system contains two very closely interacting spiral gala xies. It is k inematically very 
complex and exhibits evidence of widespread radiative shocks, which are analyzed in detail in Ric h et al.l ( l201 lb . The deprojected 
radii are plotted using the hyperleda values for the western, less obscured galaxy IC 1623 A. The eastern system is very intensely 
star-forming as seen in the infrared (e.g. Howell et al. (2010)), but is so enshrouded in dust that our optical spectra do not trace 
star-forming gas. 

IRAS F10257-4339 (NGC 3256) This advanced merger is quite nearby and well studied. As with IRAS F01053-1746, one of 
the galaxies in IRAS F10257-4339 is very extinguished. The second system and its nucl eus are revealed at longer wavelengths, 
south of the main optical nucleus ( Alonso-Herrero et al.ll2002t iRothberg & Fischerll2010l). Deprojected radii are calculated using 
the main optical nucleus. The gas and tidal tails in this system extend quite far (e.g. iRothberg & Fischerll2"010l) our IFU mosaic 
covers only the central 6 kpc, though this appears to be the physical extent of most of the ongoing star formation in this system. 
NGC 3256 also shows evidence for widespread shocks in our data and is discussed alongside IC 1623 in Rich et al. (201 1). 

IRAS F 18093-5744 N (IC 4687) We classify this galaxy as a close merger; in fact it is a member of a triplet. IC 4687 is 
undergoing a close me rger with the less massive starburst IC 4686, classed as a Wolf-Rayet galaxy bv iKovo & Continl 0999); 
lFernandeset"an (120041) . IC 4687 itself is a is also in a wide merger with the equally massive spiral IC 4689. The archived Hubble 
Space Telescope images of IC 4687 show a complex morphology tangled up with IC 4686: gas and dust from IC 4687 appear 
to be obscuring the less massive system. Our IFU data covers the entirety of IC 4686/4687 and the metallicities we measure are 
consistent with the expected metallicties in the outskirts of IC 4687 as extrapolated from the gradient we present in this paper 
as well as a low-metallicity, flattened gradient in IC 4686. The kinematic information from our data also indicates that we are 
indeed seeing gas from both systems. 

IRAS Fl 8093-5744 S (IC 4689) This s piral galaxy is sligh tly less massive and luminous than IC 4687 and less well-observed, 
though it is still intensely star forming (iHowell et al.ll2010l) . It is less morphologically disturbed than the other two interacting 
galaxies IC 4686/4687, though its gradient is quite flattened already according to our observations. Although it is widely separated 
from IC 4687, we include it as part of th e closely interac ting system of IRAS Fl 8093-5744. 

IRAS F 191 15-2124 (ESO 593-IG 008) IVaisanen et alJ (l2008i) analyze the morphology of this close pair in detail using Adaptive 
Optics K-band images. Because of the complex structure we adopt the projected radii for gradient calculations and make no 
attempt to separate the pair into two separate gradients. Given the flat values for metallicity seen, however, adapting an inclination 
and/or two disks would only lead to an even shallower calculated slope in the case of this system. 

Coalesced systems 

Our metallicity sample has two systems which are classified as coalesced mergers (cde). The most advanced mergers have 
complex orientations and morphologies which present a difficulty when calculating deprojected radii. As a conservative estimate 
we simply use the projected radii in our gradient calculations. 

IRAS 08355-4944 While Hubble Space Telescope I-band images show remnant tidal tails extending nearly 20 kpc, only the 
central 5 kpc or so appears to harbor the intense ongoing star formation in this system. Our spectra are dominated by HII 
regions, with evidence in some regions of a blue-shifted component with low-velocity shock-dominated line ratios which could 
be associated with a galactic wind. 

IRAS F10038-3338 (ESO 374-IG032) This post-merger exhibits significant ongoing star formation in its southwestern tidal arm 
unlike the two other coalesced systems in our sample. The total line emission in this region is much weaker than in the nucleus, 
inducing a large uncertainty in the extinction map and [Oil] lines, creating the discrepant values seen in the metallicity gradient 
plots for this galaxy. Either gradient, however, is consistent with the overall flattening trend seen in all of our systems. This 
system hosts an OH megamaser and has soft x-ra y emission, all consist ent wi th the advanced stage of me rging and increasingl y 
intense nuclear starburst (iHenkel & Wilsonlll990l: IStaveley-Smith et al.l ll992!: 'Darl ing & GiovaneUill200l (iwasawa et aDl2009l) . 
Our spectra also show evidence for extended off-nuclear shock er nission do minati ng in areas where there is little to no evidence 
of ongoing star formation, consistent with the IFU observations of lMonreal-Ibero et al.l (120101) . 
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Figure Al. Maps and metallicity gradients for IRAS F17222-5953. The top left two panels show [N II]/Hq& [S II]/Haline ratio maps of the system with 
Ha contours measured from our data overlaid. The top right panels show an image of the galaxy compared to the metallicity map (derived from the PP04 03N2 
diagnostic, with values scaled to the KD02 diagnostic. The black bar in the system image corresponds to a distance of one kpc. The bottom two panels show the 
run of abundance with radius calculated from the PP04 (scaled to KD()2) and the KD02 diagnostics. The least-squares fit metallicity gradient is overplotted in 
each case. The average abundance gradients for the isolated (purple) and interacting (red) sample from Rupke et al. (2010b) are plotted in the right hand panels 
for comparison. As described in the text, only points that pass our Hll-region cut in the diagnostic diagrams are plotted in the metallicity map and abundance 
gradients. 
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Figure A2. Same as fig Al for IRAS F18341-5732 




Figure A3. Same as fig Al for IRAS F13373+0105 W 
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Figure A4. Same as fig Al for IRAS F13373+0105 E 
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Figure A8. Same as fig Al for IRAS F18093-5744 S 
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Figure A9. Same as fig Al for IRAS F18093-5744 C. As discussed in the text, this appears to be a mix of IRAS F18093-5744 N and a WR galaxy, radii are 
plotted from the center of IRASF 18093-5744 and the derived values are not included in our final analysis. 
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Figure AlO. Same as fig A 1 for IRAS F191 15-2124 




Figure All. Same as fig Al for IRAS 08355-4944 
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Figure A12. Same as fig Al for IRAS F10038-3338 
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